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The ratio of B + and B° meson lifetimes was measured using data collected in 2002-2004 by the 
D0 experiment in Run II of the Fermilab Tevatron Collider. These mesons were reconstructed 
in B — > /i + vD*~ X decays, which are dominated by B°, and B — > fi + uD°X decays, which are 
dominated by B + . The ratio of lifetimes is measured to be r + /r° = 1.080±0.016 (stat)±0.014 (syst). 

PACS numbers: 14.40.Nd 



In the last few years, significant progress has been 
made, on both experimental and theoretical fronts, in 
the understanding of the lifetimes of hadrons contain- 
ing heavy quarks. Charm and bottom meson (except 
B c ) lifetimes have been measured with precisions rang- 
ing from 0.5% to 4%, although lifetimes of heavy baryons 
are not known as well |]J. On the theoretical front, pre- 
dictions are being made using a rigorous approach based 



on the heavy quark expansion (in negative powers of the 
heavy quark mass) 0, where the large mass of the bot- 
tom quark considerably simplifies calculations. Theoret- 
ical uncertainties are further reduced for ratios of life- 
times. For instance, the ratio of the B + and B° lifetimes 
has been predicted to be 1.06 ±0.02 Experimentally, 
ratios of lifetimes have smaller uncertainties, since many 
common sources of systematics cancel. 
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In this Letter, we present a measurement of the ra- 
tio of B + and B° lifetimes using a large sample of 
semileptonic B decays collected by the D0 experiment 
at Fermilab in pp collisions at y/s = 1.96 TeV. The 
data correspond to approximately 440 pb _1 of integrated 
luminosity. B mesons were selected via their decays 
B —>■ fi + vD°X Q and were classified into two exclu- 
sive groups: a "D*~" sample, containing all events with 
reconstructed D*~ — > D Q n~ decays, and a "D°" sample, 
containing all remaining events. Both simulation and 
available experimental results show that the D*~ sample 
is dominated by B° — > [i + vD*~ X decays, while the D° 
sample is dominated by B + — ► /.i + vD°X decays. 

The classification into these two samples was based on 
the presence of a slow pion from D*~ — ► D°tt~ decay, and 
thus was independent of the B-meson lifetime. Therefore, 
the ratio of the number of events in the two samples, ex- 
pressed as a function of the proper decay length, depends 
mainly on the lifetime difference between the B + and B° 
mesons. The influences of the selection criteria, detector 
properties, and some systematic uncertainties are signif- 
icantly reduced. 

The D0 detector is described in detail elsewhere 
The detector components most important to this anal- 
ysis are the central tracking and muon systems. The 
tracking system consists of a silicon microstrip tracker 
and a central fiber tracker, both located within a 2 T su- 
perconducting solenoidal magnet. The resolution for the 
distance of closest approach as provided by the tracking 
system is ~ 50 /im for tracks with pt ~ 1 GeV/c, im- 
proving asymptotically to 15 /im for tracks with px > 10 
GeV/c, where px is the component of the momentum 
perpendicular to the beam pipe. The muon system is lo- 
cated outside the calorimeters and consists of a layer of 
drift chambers and scintillation trigger counters in front 
of 1.8 T toroids, followed by two more similar layers after 
the toroids. 

Events with semi-muonic 6-hadron decays were se- 
lected using a suite of inclusive single-muon triggers in a 
three-level trigger system. Muons were identified by ex- 
trapolating tracks found in the central tracking system 
and matching them with muon track segments formed 
from hits in the muon system. Muons were required to 
have a transverse momentum p T > 2 GeV/c and total 
momentum p^ > 3 GeV/c. 

The primary vertex of the pp interaction was deter- 
mined for each event. The average position of the beam- 
collision point was included as a constraint. The preci- 
sion of the primary vertex reconstruction was on average 
about 20 /im in the plane perpendicular to the beam di- 
rection and about 40 ^m along the beam direction. 

D° candidates were selected using their D° — > K + tt~ 
decay mode. All charged particles in an event were clus- 
tered into jets using the DURHAM clustering algorithm 
with a jet pr cut-off parameter of 15 GeV/c Q. The 
D° candidate was constructed from two particles of oppo- 
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FIG. 1: (a) Invariant mass of the Kn system. The curve 
shows the result of the fit of the K + ir~ mass distribution 
with the sum of a signal Gaussian function and polynomial 
background function, (b) Mass difference Am = m(D°-ir) — 
m(D°). 



site charge belonging to the same jet as the reconstructed 
muon. Both particles were required to have pr > 0.7 
GeV/c and to form a common D° vertex. The pr of the 
D° was required to exceed 5 GeV/c. To reduce combi- 
natorial background, we required the D° vertex to have 
a positive displacement in the xy plane, relative to the 
primary vertex, with at least 4cr significance. Although 
this last requirement can bias the lifetime distribution of 
a B candidate, our analysis procedure of determining the 
ratio of B + and B° events in bins of proper time should 
remove this bias in the final result. The trajectory of the 
muon and D° candidates were required to originate from 
a common B vertex. The system was required to 

have an invariant mass between 2.3 and 5.2 GeV/c 2 . 

The masses of the kaon and pion were assigned to the 
two tracks according to the charge of the muon, assuming 
the [i + K + ir~ combination. The mass spectrum of the 
Ktt system after these selections is shown in Fig. ^a) . 
The signal in the D° peak contains 126073 ± 610 events. 

All reconstructed fj, + D° events were classified into 
three non-overlapping samples. For each /i + _D° candi- 
date, a search was made for an additional pion. The mass 
difference Am = m(D°ir) — m(D a ) for all such pions, 
when 1.8 < m(D°) < 1.9 GeV/c 2 , is shown in Fig.^b). 
The peak in this figure corresponds to the production of 
the ji + D*~ system. All events containing a pion with a 
charge opposite to that of the muon (right-charge com- 
bination) and 0.1425 < Am < 0.1490 GeV/c 2 were in- 
cluded in the D*~ R sample. All events containing a pion 
with the same charge as the muon (wrong-charge com- 
bination) and 0.1425 < Am < 0.1490 GeV/c 2 were in- 
cluded in the auxiliary D*~ w sample. This sample con- 
tains true D° but fake D*~ events and gives an estimate 
of the combinatorial background for selected n + D*~ can- 
didates. The Am distribution for such events is shown 
in Fig. ^b) as the filled histogram. All remaining events 
were assigned to the D° sample. 

Since the final (semileptonic) state has missing 
particles, including the neutrino, the proper de- 
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cay length was not determined. Instead, for 
each reconstructed candidate, the measured visible 
proper decay length x M was computed as x M — 
m B c (L T ■ p T (n+D )) /\p T (n+D°)\ 2 . L T is the vector 
in the axial plane from the primary to the i?-meson de- 
cay vertex, pt{[J. + D ) is the transverse momentum of the 
/j, + D° system and ms is the mass of the B meson, for 
which the value 5.279 GeV/c 2 was used 0. The pion 
from the D*~ decay was not used for the computation of 
the transverse momentum and the decay length. 

Candidates in each of the samples were divided into 
eight groups according to their x M value. The number of 
fi+D° events N* R (from the D*- R sample), N* w (from 
the D*~ w sample), and N® (from the D° sample) in each 
interval i (where i ranges from one to eight) were deter- 
mined from the fit of the Kir mass spectrum between 
1.72 and 2.16 GeV/c 2 with the sum of a Gaussian sig- 
nal function and a polynomial background function. The 
mean and width of the Gaussian function were fixed to 
the values obtained from the fit of the overall mass dis- 
tribution in each sample. The fitting procedure was the 
same for all samples. Table[I]gives the numbers obtained 
for each x M interval. 

The number of /j, + D*~ events for each interval i of 
x M was denned as N^i+D*-) = N* R - C ■ N* w , where 
C-N* w accounts for the combinatorial background under 
the D*~ peak as shown in Fig. ^b). The coefficient 
C = 1.27±0.03 reflects the difference in the combinatorial 
background between /j, + D°ir~ and /j, + D°tt + events. It 
was determined from the ratio of the numbers of these 
events in the interval 0.153 < Am < 0.160 GeV/c 2 . The 
number of /i + D° events in each interval i in x M was 
defined as N^D ) = N? + N* w + C ■ N* w . 

The experimental observable Vi is the ratio of /i + D*~ 
and fi + D° events in interval i of x M , i.e., r, = 
N l (fi + D*-)/N % (fi+D a ). Values of r t and statistical un- 
certainties are given in Tabic [I] The measurement of 
the lifetime difference between B + and B° is given by 
k = t + /t° — 1. It was determined from the minimization 
of x 2 (e 7r ,fc): 



X 2 (e w ,k) = J2 



(r t ~rf(^,fc)) 2 



(1) 



where rf(e n , k) is the expected ratio of n + D*~ and n + D° 
events, and is the efficiency to reconstruct the slow 
pion in the D*~ — * D Q n~ decay. e n was assumed to 
be independent of x and, along with k, was a free pa- 
rameter in the minimization. We present evidence for the 
validity of this assumption in the discussion of systematic 
uncertainties. The sum ^\ was taken over all intervals 
with positive x M . 

Information used to determine the expected ratio, 
r| (e^r , k) , included both experimental measurements as 
well as results from Monte Carlo simulations. For the 
jth B-meson decay channel, the distribution of the visible 



proper decay length (x) is given by Pj (x) — f dK Dj (K) ■ 
0(x) ■ ^-exp(— -f^)- Tj is the lifetime of the B meson, 

the if-factor, K = p!^ 5 /f>y, reflects the difference be- 
tween the observed and true momentum of the B meson, 
and 9(x) is the step function. The function Dj(K) is the 
normalized distribution of the if -factor for the jth decay 
channel. 

Transformation from the true value of x to the ex- 
perimentally measured value x M is given by fj(x M ) — 
J dx Rj(x — x M ) ■ £j(x) ■ Pj(x), where Rj(x — x M ) is the 
detector resolution, and £j(x) is the reconstruction effi- 
ciency of fj, + D° for the jth decay. It does not include 
e n for channels with D*~. Finally, the expected value 
rf(s„,k) is given by: 



F*{k) 



F?(k) + (l-e n )-F?(ky 



(2) 



Here F* ,a = J. dx M Brj ■ fj(x M ) with the summation 
J2j taken over all decays to D*~ (D°) for F* (Ff). 

For the computation of r| , the world average of the B + 
lifetime [l| was used. The B° lifetime r° was expressed 
as t° = r + /(l + k). The branching fractions B — > n + vD 
and B — * [i + vD* were taken from Ref. 0. The following 
branching fractions were derived from experimental mea- 
surements 

HSHGii : Br(B+ -> v+vD** ) = (2.67 ± 
0.37)%, Br(B+ -> n+vD**° -> D*~X) = (1.07±0.25)%, 
and Br{B° s -> n+vD* s *-) = (2.3t 2 ;3)%. D ** states in- 
clude both narrow and wide D** resonances and non- 
resonant DX and D*X production. Regarding the pos- 
sible decays of D**~ , there is no experimental data on 
the Br(D**~ — * D*~X). Its central value was therefore 
set to 0.35 and it was varied between 0.0 and 1.0 to es- 
timate the systematic uncertainty from this source. All 
other branching fractions were derived assuming isospin 
invariance. 

The distributions Dj(K), Rj(x), and £j(x) were taken 
from the Monte Carlo simulation, and the corresponding 
systematic uncertainties were taken into account. All 
processes involving b hadrons were simulated using the 
EvtGen 0] generator interfaced to pythia Q and fol- 
lowed by the full modeling of the detector response and 
event reconstruction. The semileptonic 6-hadron decays 
were generated using the ISGW2 model [l^ . 

Assuming the given branching fractions and recon- 
struction efficiencies, the decay B — > fi + D*~X contains 
(89 ±3)% B°, (10 ±3)% B+, and (1±1)% B° s , while the 
decay B -> ^+D°X contains (83 ± 3)% B+, (15 ± 4)% 
5°, and (2± 1)% B° s . 

A special study showed that in addition to the main de- 
cay B -► Ijl+D q X, the decay B -► t+D°X -> n + vvD Q X 
results in a (5 ± 2)% contribution and the process cc — > 
fi+D°X a (10 ± 7)% contribution to the selected fi + D° 
sample. These processes were taken into account in the 
analysis. 
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TABLE I: Definition of the intervals in visible proper decay length, x . For each interval i, the number of events in the D* 
D*~ w and D° samples, the ratio n, and the expected value r\ for r + /r° — 1 = 0.080 are given. 



i 


x M range (cm) 


N* R 


N* w 


N? 






1 


-0.1 - 0.0 


1714 ± 53 


89 ± 22 


5225 ± 151 


0.295 ± 0.015 


0.309 


2 


0.0 - 0.02 


6213 ± 94 


200 ± 28 


18134 ± 222 


0.321 ± 0.007 


0.315 


3 


0.02 - 0.04 


5941 ± 91 


169 ± 22 


17703 ± 208 


0.317 ± 0.007 


0.313 


4 


0.04 - 0.07 


6424 ± 94 


213 ± 23 


19707 ± 216 


0.305 ± 0.006 


0.308 


5 


0.07 - 0.10 


4029 ± 74 


115 ± 17 


12885 ± 171 


0.295 ± 0.007 


0.300 


6 


0.10 - 0.15 


3459 ± 68 


106 ± 16 


11532 ± 162 


0.282 ± 0.007 


0.291 


7 


0.15 - 0.25 


2253 ± 57 


58 ± 13 


7567 ± 137 


0.283 ± 0.009 


0.276 


8 


0.25 - 0.40 


518 ± 28 


2 ± 6 


1875 ± 75 


0.274 ± 0.019 


0.256 




0.25 



-0.1 0.1 0.2 0.3 0.4 

Visible Proper Decay Length (cm) 

FIG. 2: Points with the error bars show the ratio of the num- 
ber of events in the [i + D*~ and jj, + D° samples as a function 
of the visible proper decay length. The result of the mini- 
mization of Eq. with k = 0.080 is shown as a histogram. 



Using all these inputs, the minimization of the x 2 dis- 
tribution, Eq. 0, gives: k = t+/t° - 1 = 0.080 ± 
0.016 (stat). The \ 2 & t the minimum is 4.2 for 5 d.o.f, 
e„ is 0.864 ± 0.006 (stat), and the global correlation co- 
efficient between k and e n is 0.18. The simulation pre- 
dicted En = 0.877 ± 0.003. The reasonable agreement in 
Ett between data and simulation reflects good consistency 
of input efficiencies and branching fractions with exper- 
imental data. Figure [21 presents the T{ values together 
with the result of the fit. 

The influence of various sources of systematic uncer- 
tainty on the final result is summarized in Table Irfl Dif- 
ferent contributions can be divided into three groups. 
The first part includes uncertainties coming from the ex- 
perimental measurements, e.g., branching fractions and 
lifetimes. All inputs were varied by one standard devia- 
tion. Only the most significant contributions are listed as 
individual entries in Table [n] all remaining uncertainties 
are combined into a single entry "other contributions." 

The second group includes uncertainties due to the in- 
puts taken from the Monte Carlo simulation. They were 
estimated as follows. The uncertainty due to the decay 
length dependence of the efficiencies e(B — > fi + vD a X) 
was obtained by repeating the analysis with decay length 



independent efficiencies used for all decay modes. This 
dependence almost cancels in the ratio of the number of 
events in the two samples, leading to the reduced sys- 
tematic uncertainty in r + /r°. 

The variation of the efficiency from channel to channel 
arises from differences in the kinematics of S-meson de- 
cays and thus depends on their modeling in simulation. 
To estimate the uncertainty in the efficiency due to this 
effect, an alternative HQET model of B- meson decays 
|l3| was implemented, and the selection cuts on the pt 
of the /i + and D° were varied over a wide range. 

The same alternative model and the variation of pt 
cuts were used to study the model dependence of the K- 
factors. In all cases, the variation of the average value 
of i-T-factors did not exceed 2%. Distributions of K- 
factors were determined separately for B — > [jl + vD®, 
B -> n+vD*, B -> n+vD** -> D°X, and B -> 
H + vD** — ► D*X. To estimate the uncertainty due to 
the modeling of D** decays, which include both reso- 
nant and non-resonant components and are not yet well 
understood, the analysis was repeated with the distribu- 
tions of if-factors from B — > D** — > D°(D*) decays set 
to be the same as for B — > D°(D*) decays. 

The selection of the slow pion was made indepen- 
dently of the B lifetime, and the efficiency e w was as- 
sumed constant in the minimization. A dedicated study 
of Ks — > it + it~ decays showed good stability of the track 
reconstruction efficiency with the change of decay length 
over a wide range. The slope in the efficiency was esti- 
mated to be 0.0038 ± 0.0059 cm" 1 . The independence 
of on the decay length was also verified in simulation, 
where no deviation from the constant value was detected 
within available simulation statistics. 

The average decay length resolution, approximately 35 
/im for this measurement, and the fraction of events with 
larger resolution, modeled by a Gaussian function with 
resolution of 1700 ^m, were varied over a wide range, sig- 
nificantly exceeding the estimated difference in resolution 
between data and simulation. 

The ratio of events with negative decay length in D*~ 
and D° samples (the first row in Table is sensitive to 
the differences in resolution of these two samples. The 
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TABLE II: Summary of systematic uncertainties. 



Source 


A(r+/r°) 


Br(B IJ -» fi+uD*-) 


0.0005 


Br(B + -> [i + vD* ) 


0.0010 


Br(B + -» (i + vD u ) 


0.0009 


Br(B + -» (i + vD*-tt + X) 


0.0059 


Br(B° -> fi + vD~X) 


0.0009 


DT'^D'-X 


0.0020 


cc — > fi + uD°X contribution 


0.0015 


Other contributions 


0.0006 


e(B — > fi + uD u X), decay length dependence 


0.0014 


e(_B — » (i + vD° X), average value 


0.0030 


e,r, decay length dependence 


0.0036 


decay length resolution 


0.0024 


Difference in D*~ and _D° resolution 


0.0053 


AT-factors, average value 


0.0032 


AT-factors, difference between channels 


0.0013 


Fitting procedure 


0.0086 


Background level under D*~ 


0.0004 


Total 


0.0136 



comparison of this ratio in data with the simulation was 
used to estimate the systematic uncertainty due to this 
difference. 

Since the /j, + D*~ and fj, + D° event yields as a function 
of proper decay length are extracted by fitting the re- 
spective mass distributions with signal and background 
functional forms, the fitting procedure can be another 
source of systematic uncertainty. Different parameter- 
izations of signal and background functions were used. 
The maximal variation of the result obtained was taken 
as the systematic uncertainty due to this source. Finally, 
the uncertainty in the background level under the D*~ 
peak in Fig.QIb) was also taken into account. 

Estimated systematic uncertainties from different 
sources were added in quadrature and the total system- 
atic uncertainty on the ratio of lifetimes is <t(t + /t°) 
0.0 1 4 Various consistency checks of this measurement 
were also performed. The total sample of events was di- 
vided into two parts using different criteria, such as the 
sign of the muon rapidity, polarity of the solenoid, charge 
of the muon, pt of the muon, position of the primary in- 
teraction, etc. The measurement was repeated indepen- 
dently for each sample. The definition of proper decay 
length intervals was varied, one more interval, 0.4 - 0.8 
cm, was added, and the last interval, 0.25 - 0.4 cm, was 
removed from the fit. In all cases, the results are con- 
sistent within statistical uncertainties. Finally, the mea- 
surement of the ratio of lifetimes was performed with sim- 
ulated events. The resulting value fc MC = 0.084 ± 0.015 
agrees well with the generated lifetime ratio fc gcn = 0.070. 

In summary, the ratio of B + and B meson lifetimes 
was found to be: 

r+ 

k = — - 1 = 0.080 ± 0.016 (stat) ± 0.014 (syst). (3) 



This result is the most precise measurement of this pa- 
rameter, and agrees well with the world average value 
k = 0.086 ± 0.017 0. Improved precision of the ratio of 
B + and B° lifetimes will allow a better test of theoreti- 
cal predictions, especially those inputs to the calculations 
that rely on lattice QCD or on other non-perturbative 
methods HQ. 
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